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o ■ 

The status of the constrained minimal supersymmetric standard model 
O ' (CMSSM) will be discussed in light of our current understanding of the relic 

. density after WMAP. A global likelihood analysis of the model is performed 

including data from direct Higgs searches, global fits to electroweak data, 
CO . b sj, the anomalous magnetic moment of the muon, as well as the cos- 

mological relic density. Also considered are models which relax and further 
. , constrain the CMSSM. Prospects for dark matter detection in colliders and 

' cryogenic detectors will be briefly discussed. 

o , 

^ . 1 Introduction 
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■ Supersymmetric models with conserved i?-parity contain one new stable par- 

ticle which is a candidate for cold dark matter (CDM) 1 . There are very 
strong constraints, however, forbidding the existence of stable or long lived 
particles which are not color and electrically neutral. The sneutrino 2 is one 
possible candidate, but in the MSSM, it has been excluded as a dark matter 
(— I I candidate by direct T and indirect W- searches. Another possibility is the 

gravitino and is probably the most difficult to exclude. This possibility has 
been discussed recently in the CMSSM context 0. I will concentrate on the 
remaining possibility in the MSSM, namely the neutralinos. 
5— i ' There are four neutralinos, each of which is a linear combination of the 

^ i? = — 1, neutral fermions the wino W'^, the partner of the 3rd component 

of the SU{2)l gauge boson; the bino, B, the partner of the U{1)y gauge 
boson; and the two neutral Higgsinos, Hi and H2. In general, the neutralino 
mass eigenstates can be expressed as a linear combination 

X^aB + (3W^ + -yHi + 6H2 (1) 

The solution for the coefficients a,/?, 7 and S for neutralinos that make up 
the LSP can be found by diagonalizing the mass matrix which depends on 
Mi{M2) which are the soft supersymmetry breaking U(l) (SU(2)) gaugino 
mass terms, fj,, the supersymmetric Higgs mixing mass parameter and the 
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two Higgs vacuum expectation values, vi and V2- One combination of these is 
related to the Z mass, and therefore is not a free parameter, while the other 
combination, the ratio of the two vevs, tan P, is free. 

The most general version of the MSSM, despite its minimality in parti- 
cles and interactions contains well over a hundred new parameters. The study 
of such a model would be untenable were it not for some (well motivated) 
assumptions. These have to do with the parameters associated with super- 
symmetry breaking. It is often assumed that, at some unification scale, all 
of the gaugino masses receive a common mass, mi/2- The gaugino masses 
at the weak scale are determined by running a set of renormalization group 
equations. Similarly, one often assumes that all scalars receive a common 
mass, mo, at the GUT scale. These too are run down to the weak scale. The 
remaining supersymmetry breaking parameters are the trilinear mass terms, 
Aq, which I will also assume are unified at the GUT scale, and the bilinear 
mass term B. There are, in addition, two physical CP violating phases which 
will not be considered here. 

The natural boundary conditions at the GUT scale for the MSSM would 
include fj, and B in addition to mi/2, mo, and Ao- In this case, upon running 
the RGEs down to a low energy scale and minimizing the Higgs potential, 
one would predict the values of Mz, tan/3 (in addition to all of the sparticle 
masses). Since Mz is known, it is more useful to analyze supersymmetric 
models where Mz is input rather than output. It is also common to treat 
tan/3 as an input parameter. This can be done at the expense of shifting ^ 
(up to a sign) and B from inputs to outputs. This model is often referred 
to as the constrained MSSM or CMSSM. Once these parameters are set, the 
entire spectrum of sparticle masses at the weak scale can be calculated. In 
the CMSSM, the solutions for /i generally lead to a neutralino which which 
very nearly a pure B. 

2 The CMSSM after WMAP 

For a given value of tan /3, Aq, and sgn{ii), the resulting regions of acceptable 
relic density and which satisfy the phenomcnological constraints can be dis- 
played on the mi/2 — Too plane. In Fig.ni, the light shaded region corresponds 
to that portion of the CMSSM plane with tan/3 = 10, Ao = 0, and /i > 
such that the computed relic density yields 0.1 < J7^/i^ < 0.3. At relatively 
low values of mi/2 ^^d toq, there is a large 'bulk' region which tapers off 
as my 2 is increased. At higher values of mg, annihilation cross sections are 
too small to maintain an acceptable relic density and Q^h^ > 0.3. Although 
sfermion masses are also enhanced at large my2 (due to RGB running), co- 
annihilation processes between the LSP and the next lightest sparticle (in 
this case the Ti) enhance the annihilation cross section and reduce the relic 
density. This occurs when the LSP and NLSP are nearly degenerate in mass. 
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The dark shaded region has m~ < and is excluded. Neglecting coannihi- 
lations, one would find an upper bound of 450 GeV on my2i corresponding 
to an upper bound of roughly 200 GeV on . The effect of coannihilations is 
to create an allowed band about 25-50 GeV wide in mo for TOi/2 ^ 1400 GeV, 
which tracks above the mf^ = contour 




Fig. 1. The (mi/2,n^o) planes for (a) tan/3 = 10 and /j, > 0, assuming Ao = 

0,mt — 175 GeV and mb{mb)s'M ~ 4-25 GeV. The near-vertical (red) dot-dashed 
lines are the contours mu = 114 GeV, and the near-vertical (black) dashed line 
is the contour rn^± = 104 GeV. Also shown by the dot-dashed curve in the lower 
left is the comer excluded by the LEP bound of > 99 GeV. The medium (dark 
green) shaded region is excluded by b —> 57, and the light (turquoise) shaded area 
is the cosmologically preferred regions with 0.1 < il^h^ < 0.3. In the dark (brick 
red) shaded region, the LSP is the charged fi. The region allowed by the E821 
measurement of at the 2-a level, is shaded (pink) and bounded by solid black lines, 
with dashed lines indicating the 1-a ranges. In (b), the relic density is restricted to 
the range 0.094 < Q-^h^ < 0.129. 



Also shown in Fig. ^ are the relevant phcnomcnological constraints. 
These include the limit on the chargino mass: rn^± > 104 GeV [7], on the 
selectron mass: mg > 99 GeV ||8| and on the Higgs mass: ruh > 114 GeV 
The former two constrain mi/2 and mo directly via the sparticle masses, and 
the latter indirectly via the sensitivity of radiative corrections to the Higgs 
mass to the sparticle masses, principally m^i. FeynHiggs ^U] is used for the 
calculation of m/i. The Higgs limit imposes important constraints principally 
on mi/2 particularly at low tan/3. Another constraint is the requirement that 
the branching ratio for & — > 57 is consistent with the experimental mea- 
surements 22 . These measurements agree with the Standard Model, and 
therefore provide bounds on MSSM particles |12[ I13| . such as the chargino 
and charged Higgs masses, in particular. Typically, the & — *■ 57 constraint 
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is more important for /i < 0, but it is also relevant for /i > 0, particularly 
when tan/3 is large. The constraint imposed by measurements of 6 — > S7 also 
excludes small values of Finally, there are regions of the (toi/2,"2o) 

plane that are favoured by the BNL measurement of 5^ — 2 at the 2-cr 
level, corresponding to a deviation from the Standard Model calculation ^S] 
using e^e~ data. One should be however aware that this constraint is still 
under active discussion. 

The preferred range of the relic LSP density has been altered significantly 
by the recent improved determination of the allowable range of the cold 
dark matter density obtained by combining WMAP and other cosmological 
data: 0.094 < Qcdm < 0.129 at the 2-cr level [TBj. In the second panel of 
Fig. nj we see the effect of imposing the WMAP range on the neutralino 
density ^\ • We see immediately that (i) the cosmological regions are 

generally much narrower, and (ii) the 'bulk' regions at small mi/2 and ttiq 
have almost disappeared, in particular when the laboratory constraints are 
imposed. Looking more closely at the coannihilation regions, we see that (iii) 
they are significantly truncated as well as becoming much narrower, since 
the reduced upper bound on flyji^ moves the tip where = rrif to smaller 
mi/2 so that the upper limit is now mi/2 ^ 950 GeV or < 400 GeV. 



tanB = 50, n>0 tan B = 50 , Li > 

L L-L 1 linn L LJI 
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Fig. 2. As m Fig.[Zfor tan/3 = 50. 

Another mechanism for extending the allowed CMSSM region to large 
is rapid annihilation via a direct-channel pole when ~ ^ruA [^Oll^- Since 
the heavy scalar and pseudoscalar Higgs masses decrease as tan /3 increases, 
eventually 2m^ ~ ttia yielding a 'funnel' extending to large mi/2 and mo 
at large tan /?, as seen in the high tan (3 strips of Fig. [3 As one can see, 
the impact of the Higgs mass constraint is reduced (relative to the case with 
tan f3 = 10) while that of 6 ^ S7 is enhanced. 
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Shown in Fig.Oare the WMAP hnes 17 of the {mi/2,mo) plane allowed 
by the new cosmological constraint 0.094 < H^h^ < 0.129 and the laboratory 
constraints listed above, for /i > and values of tan/3 from 5 to 55, in 
steps Z\(tan/3) = 5. We notice immediately that the strips are considerably 
narrower than the spacing between them, though any intermediate point in 
the (mi/2,'7^o) plane would be compatible with some intermediate value of 
tan p. The right (left) ends of the strips correspond to the maximal (minimal) 
allowed values of mi/2 and hence m^. The lower bounds on mi/2 are due to 
the Higgs mass constraint for tan/3 < 23, but are determined by the & — > 57 
constraint for higher values of tan/3. 
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Fig. 3. The strips display the regions of the (mi/2,mo) plane that are compatible 
with 0.094 < fi-^h^ < 0.129 and the laboratory constraints for fi > and tan/3 = 
5, 10, 15, 20, 25, 30, 35, 40, 45, 50, 55. The parts of the strips compatible with 3^ — 2 
at the 2-a level have darker shading. 



Finally, there is one additional region of acceptable relic density known 
as the focus-point region |22| . which is found at very high values of uiq. An 
example showing this region is found in Fig.^ plotted for tan/3 = 10, /i > 0, 
and mt = 175 TeV. As mo is increased, the solution for fi at low energies 
as determined by the electroweak symmetry breaking conditions eventually 
begins to drop. When ^ < mi/2, the composition of the LSP gains a strong 
Higgsino component and as such the relic density begins to drop precipitously. 
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These effects are both shown in Fig. |51 where the value of /i and f2h^ are 
plotted as a function of mo for fixed mi/2 = 300 GeV and tan (3 = 10. As mp 
is increased further, there are no longer any solutions for fi. This occurs in 
the shaded region in the upper left corner of Fig. 0] 



tan B = 10 , |i > 

5000 I 




100 200 300 400 500 600 700 800 900 1000 



mi/2 (GeV) 

Fig. 4. As in Fig.^p., where the range in mo is extended to 5 TeV. In the shaded 
region at very high mo, there are no solutions for which respect the low energy 
electroweak symmetry breaking conditions. 

Fig. 13 also exemplifies the degree of fine timing associated with the focus- 
point region. While the position of the focus-point region in the mo,TOi/2 
plane is not overly sensitive to supersymmetric parameters, it is highly sensi- 
tive to the top quark Yukawa coupling which contributes to the evolution of 
fi 1^3121]. As one can see in the figure, a change in mt of 3 GeV produces a 
shift of about 2.5 TeV in tuq. Note that the position of the focus-point region 
is also highly sensitive to the value of Aq/itiq. In Fig.jS] Aq — was chosen. 
For Ao/niQ = 0.5, the focus point shifts from 2.5 to 4.5 TeV and moves to 
larger mg as Ao/mo is increased. 

3 A Likelihood analysis of the CMSSM 

Up to now, in displaying acceptable regions of cosmological density in the 
mo, nil /2 plane, it has been assumed that the input parameters are known 
with perfect accuracy so that the relic density can be calculated precisely. 
While all of the beyond the standard model parameters are completely un- 
known and therefore carry no formal uncertainties, standard model parame- 
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Fig. 5. The value of ji as a function of mo for fixed m-^/ 2 ~ 300 GeV and tan (3 — 10 
for two choices of rut as indicated. The scale on the right gives the value of Qh^ . 
The curves corresponding to this is scale rise sharply at low mo to values much 
larger than 1. For rut = 175 GeV and mo ~ 2500 GeV, the value of fih^ drops to 
acceptable values when /i becomes small. When the mt — 178 GeV, Qh^ drops at 
mo « 5000 GeV. 

ters such as the top and bottom Yukawa couphngs are known but do carry 
significant uncertainties. Indeed, we saw that in the case of the focus-point 
region, there is an intense sensitivity of the rehc density to the top quark 
Yukawa. Other regions in the r7io,"^i/2 plane, such as those corresponding 
to the rapid annihilation funnels are also very sensitive to the 3rd generation 
Yukawas. 

The optimal way to combine the various constraints (both phenomenolog- 
ical and cosmological) is via a likelihood analysis, as has been done by some 
authors both before [23 and after ^3 the WMAP data was released. When 
performing such an analysis, in addition to the formal experimental errors, it 
is also essential to take into account theoretical errors, which introduce sys- 
tematic uncertainties that are frequently non-negligible. Recently, we have 
preformed an extensive likelihood analysis of the CMSSM [221 ■ Included is 
the full likelihood function for the LEP Higgs search, as released by the LEP 
Higgs Working Group. This includes the small enhancement in the likelihood 
just beyond the formal limit due to the LEP Higgs signal reported late in 
2000. This was re-evaluated most recently in [5], and cannot be regarded as 
significant evidence for a light Higgs boson. We have also taken into account 
the indirect information on provided by a global fit to the precision elec- 
troweak data. The likelihood function from this indirect source does not vary 
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rapidly over the range of Higgs masses found in the CMSSM, but we included 
this contribution with the aim of completeness. 

The interpretation of the combined Higgs likelihood, C^xpi in the (?tI]^/2i "^o) 
plane depends on uncertainties in the theoretical calculation of rrih. These in- 
clude the experimental error in rrit and (particularly at large tan/3) mf,, and 
theoretical uncertainties associated with higher-order corrections to mt ■ Our 
default assumptions are that mt = 175 ± 5 GeV for the pole mass, and 
mb = 4.25 ± 0.25 GeV for the running MS mass evaluated at mf, itself. 
The theoretical uncertainty in to/j, ath, is dominated by the experimental 
uncertainties in rritfi, which are treated as uncorrclatcd Gaussian errors: 

Typically, we find that {dmij dnit) ~ 0.5, so that ath is roughly 2-3 GeV. 

The combined experimental likelihood, Cexp, from direct searches at 
LEP 2 and a global electroweak fit is then convolved with a theoretical like- 
lihood (taken as a Gaussian) with uncertainty given by ath from (0) above. 
Thus, we define the total Higgs likelihood function, Ch, as 

Chimh) = / dm'^ Cexpim'^) e-^-™'')'/^-- , (3) 



where A/" is a factor that normalizes the experimental likelihood distribution. 

In addition to the Higgs likelihood function, we have included the like- 
lihood function based on 6 ^ 57. The branching ratio for these decays has 
been measured by the CLEO, BELLE and BaBar collaborations and 
we took as the combined value Bib sj) = (3.54 ± 0.41 ± 0.26) x 10""*. 
The theoretical prediction |12[ [T^ contains uncertainties which stem from 
the uncertainties in mb, a^, the measurement of the semileptonic branching 
ratio of the B meson as well as the effect of the scale dependence. While the 
likelihood function based on the measurements of the anomalous magnetic 
moment of the muon was considered in |26| . it will not be discussed here. 

Finally, in calculating the likelihood of the CDM density, we take into 
account the contribution of the uncertainties in nit^b- We will see that the 
theoretical uncertainty plays a very significant role in this analysis. The like- 
lihood for f2h^ is therefore, 

= e ^ ' ' , (4) 

v27r(T 



where = al^^ + cr^^ , with aexp taken from the WMAP [I^ result and (Tj^ 
from 101, replacing rrih by fih? . 

The total likelihood function is computed by combining all the compo- 
nents described above: 



(5) 
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The likelihood function in the CMSSM can be considered a function of two 
variables, £(ot('Tii/2, Jtiq), where and toq are the unified GUT-scale 

gaugino and scalar masses respectively. Results are based on a Bayesian anal- 
ysis, in which a prior range for mi/2 is introduced in order to normalize the 
conditional probability obtained from the likelihood function using Bayes' 
theorem. Although it is possible to motivate some upper limit on toi/2, e.g., 
on the basis of naturalness |27l 1241 |2HI , one cannot quantify any such limit 
rigorously. Within the selected range, we adopt a flat prior distribution for 
mi/2, and normalize the volume integral: 



for each value of tan /3, combining where appropriate both signs of fi. We note 
that no such prior need be specified for mo. For any given value of TOi/2, the 
theory is well defined only up to some maximum value of mo, above which 
radiative electroweak symmetry breaking is no longer possible. We always 
integrate up to that point, adopting also a flat prior distribution for mo. 

In Fig.Elthe likelihood along slices through the CMSSM parameter space 
for tan/3 = 10, = 0, ^ > 0, and mi/2 — 300 and 800 GeV is shown in 
the left and right panels, respectively, plotting the likelihood as a function 
of mg. The solid red curves show the total likelihood function calculated 
including the uncertainties which stem from the experimental errors in rut 
and mi,. The peak at low mo is due to the coannihilation region. The peak 
at mo ~ 2500(4500) GeV for mi/2 = 300(800) GeV is due to the focus-point 
region. Also shown in Fig.Elare the 68%, 90%, and 95% CL (horizontal) lines, 
corresponding to the iso-likelihood values of the fully integrated likelihood 
function corresponding to the solid (red) curve. 

The focus-point peak is suppressed relative to the coannihilation peak at 
low mo because of the theoretical sensitivity to the experimental uncertainty 
in the top mass. We recall that the likelihood function is proportional to cr~^, 
and that a which scales with d{fl^h'^) / drrit, is very large at large mo |211- The 
error due to the uncertainty in m* is far greater in the focus-point region than 
in the coannihilation region. Thus, even though the exponential in Cq^i^ is 
of order unity near the focus-point region when ~ 0.1, the prefactor is 

very small due the large uncertainty in the top mass. This accounts for the 
factor of > 1000 suppression seen in Fig. El when comparing the two peaks of 
the solid red curves. 

We note also that there is another broad, low-lying peak at intermediate 
values of mo. This is due to a combination of the effects of a in the pref- 
actor and the exponential. We expect a bump to occur when the Gaussian 
exponential is of order unity, i.e., ^ \/2Amt dQ^h? / drnt- f2^h? ~ 10 

at large mo for our nominal value m* = 175 GeV, but it varies significantly 
as one samples the favoured range of mt within its present uncertainty. The 
competition between the exponential and the prefactor would require a large 
theoretical uncertainty in fi^h?: dH^h^ /dnit ~ 2 for Amt = 5 GeV. This 
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Fig. 6. The likelihood function along slices in mo through the CMSSM parameter 
space for tan/3 = 10, Aq — 0, ^ > and mi/2 = 300, 800 GeV in the left and right 
panels, respectively. The red (solid) curves are calculated using the current errors 
in mt and rUb, the green dashed curve with no error in rut, the violet dotted lines 
with Arrit = 0.5 GeV, and the blue dashed-dotted lines with Arrit = 1 GeV. 



occurs when mo ^ 1000 GeV, which is the position of the broad secondary 
peak in Fig. At higher mo, cr continues to grow, and the prefactor sup- 
presses the likelihood function until fi^^h^ drops to ^ 0.1 in the focus-point 
region. 

As is clear from the above discussion, the impact of the present exper- 
imental error in nit is particularly important in this region. This point is 
further demonstrated by the differences between the curves in each panel, 
where we decrease ad hoc the experimental uncertainty in mt- As Arrit is 
decreased, the intermediate bump blends into the broad focus-point peak. 
When the uncertainties in nit and rrib are set to 0, we obtain a narrow peak 
in the focus-point region. 

Using the fully normalized likelihood function Ctot obtained by combining 
both signs of fi for each value of tan /3, we can determine the regions in 
the {niif2,'mo) planes which correspond to specific CLs. Fig. [3 extends the 
previous analysis to the entire (1111/2, fno) plane for tan/3 = 10 and Aq = 0, 
including both signs of /i. The darkest (blue), intermediate (red) and lightest 
(green) shaded regions are, respectively, those where the likelihood is above 
68%, above 90%, and above 95%. Overall, the likelihood for /i < is less than 
that for /i > due to the Higgs and 6 ^ 57 constraints. Only the bulk and 
coannihilation-tail regions appear above the 68% level, but the focus-point 
region appears above the 90% level, and so cannot be excluded. 

The bulk region is more apparent in the right panel of Fig. [71 for /i > 
than it would be if the experimental error in mt and the theoretical error 
in mti were neglected. Fig. |S1 complements the previous figures by showing 
the likelihood functions as they would appear if there were no uncertainty in 
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Fig. 7. Contours of the likelihood at the 68%, 90% and 95% levels for tan (5 = 10, 
Ao = and /j. > (left panel) or /i < (right panel), calculated using information 
of ruh, b ^ sj and ficDMh^ and the current uncertainties in mt and mt,. 



mt, keeping the other inputs the same. We see that, in this case, both the 
coannihilation and focus-point strips rise above the 68% CL. 
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Fig. 8. As in Fig. Y\ but assuming zero uncertainty in mt. 



Fig. El shows the Ukehhood projection for tan/3 = 50, — Q and ^ > 
0. In this case, regions at small mi/2 and toq are disfavoured by the b 
sj constraint. The coannihilation region is broadened by a merger with the 
rapid-annihilation funnel. Both the coannihilation and the focus-point regions 
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feature strips allowed at the 68% CL, and these are linked by a bridge at the 
95% CL. 
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Fig. 9. Likelihood contours as in Fig. but for tan/3 = 50, Aq = and > 0. 



4 Beyond the CMSSM 

The results of the CMSSM described in the previous sections are based heav- 
ily on the assumptions of universality of the supersymmetry breaking param- 
eters. One of the simplest generalizations of this model relaxes the assumption 
of universality of the Higgs soft masses and is known as the NUHM [2^1 In 
this case, the input parameters include fi and ttia , in addition to the standard 
CMSSM inputs. In order to switch /i and ruA from outputs to inputs, the two 
soft Higgs masses, toi, m2 can no longer be set equal to toq and instead are 
calculated from the electroweak symmetry breaking conditions. The NUHM 
parameter space was recently analyzed 20; and a sample of the results are 
shown in Fig. 1101 

In the left panel of Fig. 1101 we see a nT.x/2i"^o plane with a relative low 
value of fi. In this case, an allowed region is found when the LSP contains 
a non-negligible Higgsino component which moderates the relic density in- 
dependent of mp. To the right of this region, the relic density is too small. 
In the right panel, we see an example of the ttia, ^J■ plane. The crosses corre- 
spond to CMSSM points. In this single pane, we see examples of acceptable 
cosmological regions corresponding to the bulk region, co-annihilation region 
and s-channel annihilation through the Higgs pseudo scalar. 

Rather than relax the CMSSM, it is in fact possible to further con- 
strain the model. While the CMSSM models described above are certainly 
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my2 (GeV) |x (GeV) 



Fig. 10. a) The NUHM (mi/2,mo) plane for tan/3 = 35, (a) fi = 400 GeV and 
VIA = 700 GeV b)the NUHM {p,mA) plane for tan/3 = 10, mo = 100 GeV 
and mi/2 = 300 GeV, with Aq = 0. The (red) dot-dashed lines are the con- 
tours ruh = 114 GeV, and the near-vertical (black) dashed lines are the contours 
m^± = 103.5 GeV. The dark (black) dot-dashed lines indicate the GUT stability 
constraint. Only the areas inside these curves (small /j.) are allowed by this con- 
straint. The light (turquoise) shaded areas are the cosmologically preferred regions 
with 0.1 < Q^h? < 0.3. The darker (blue) portion of this region corresponds to the 
newer WMAP densities. The dark (brick red) shaded regions is excluded because a 
charged particle is lighter than the neutralino, and the lighter (yellow) shaded re- 
gions is excluded because the LSP is a sneutrino. The medium (green) shaded region 
is excluded by b —> 37. The regions allowed by the g — 2 constraint are shaded (pink) 
and bounded by solid black lines. The solid (blue) curves correspond to — mA/2. 

mSUGRA inspired, minimal supergravity models can be argued to be still 
more predictive. Let us assume that supersymmetry is broken in a hidden 
sector so that the superpotential can be written as a sum of two terms, 
W = F{(f>) + g{Q, where ^ represents all observable fields and C all hidden 
sector fields. We furthermore must choose g{() such that when ( picks up 
a vacuum expectation value, supersymmetry is broken. When the potential 
is expanded and terms inversely proportional to Planck mass are dropped, 
one finds |3(Jj 1) scalar mass universality with mo = (g); 2) trilinear mass 
universality with Aq = {dg/d(){Q + {g){C)^; and 3) Bq = Aq — rriQ. 

In the simplest version of the theory [SJ, the universal trilinear soft 
supersymmetry-breaking terms are ^ = (3 — •\/3)too and bilinear soft 
supersymmetry-breaking term is i? = (2 — V3)wo, i.e., a special case of 
the general relation above between B and A. 

Given a relation between Bq and we can no longer use the standard 
CMSSM boundary conditions, in which toi/2, tuq, Aq, tan/3, and sgn{^) 
are input at the GUT scale with /i and B determined by the electroweak 
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symmetry breaking condition. Now, one is forced to input Bq and instead 
tan f3 is calculated from the minimization of the Higgs potential [32] ■ 

In Fig. 111! the contours of tan f3 (solid blue hues) in the (mi/27 ruo) planes 
for two values oi A = Ao/mo, B = Bo/niQ = A — 1 and the sign of fj, are 
displayed Also shown are the contours where m^± > 104 GeV (near- 
vertical black dashed lines) and rrih > 114 GeV (diagonal red dash-dotted 
lines). The excluded regions where > nif-^ have dark (red) shading, those 
excluded by 6 ^ 57 have medium (green) shading, and those where the relic 
density of ncutrahnos lies within the WMAP range 0.094 < H^h"^ < 0.129 
have light (turquoise) shading. Finally, the regions favoured by 17^ — 2 at the 
2-tT level are medium (pink) shaded. 



A = 3-\l3;H>0 A=2;n>0 
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Fig. 11. Examples 0/ (mi/2, tiq) planes with contours of tan P superposed, for 
H > and (a) the simplest Polonyi model with A — 3 — \/3, -B — A — 1 and (h) 
A = 2.0, ]3 — A — 1. In each panel, we show the regions excluded by the LEP lower 
limits on MSSM particles, those ruled out by b —> s'y decay (medium green shading), 
and those excluded because the LSP would be charged (dark red shading). The region 
favoured by the WMAP range has light turquoise shading. The region suggested by 
~ 2 IS medium (pink) shaded. 

In panel (a) of Fig. ^2 we see that the Higgs constraint combined with 
the relic density requires tan/3 > 11, whilst the relic density also enforces 
tan/3 < 20. For a given point in the toi/2 — plane, the calculated value 
of tan (3 increases as A increases. This is seen in panel (b) of Fig. when 
A = 2.0, close to its maximal value for ^ > 0, the tan/3 contours turn over 
towards smaller mi/2, and only relatively large values 25 <, tan/3 ^ 35 are 
allowed by the b ^ and QcDMh"^ constraints, respectively. For any given 
value of A, there is only a relatively narrow range allowed for tan/3. 
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5 Detectability 

The question of detectability with respect to supersymmetric models is of 
key importance particularly with the approaching start of the LHC. As an 
aid to the assessment of the prospects for detecting sparticles at different 
accelerators, benchmark sets of supersymmetric parameters have often been 
found useful, since they provide a focus for concentrated discussion 1^ 
A set of proposed post-LEP benchmark scenarios [HSI were chosen to 
span the CMSSM. Five of the chosen points are in the 'bulk' region at small 
mi/2 and mo, four are spread along the coannihilation 'tail' at larger toi/2 
for various values of tan (3. Two points are in rapid-annihilation 'funnels' at 
large mi/2 and itiq. Two points were chosen in the focus-point region at large 
mo- The proposed points range over the allowed values of tan/3 between 5 
and 50. 

In Fig. El a comparison of the numbers of different MSSM particles that 
should be observable at different accelerators in the various benchmark sce- 
narios ISH], ordered by their consistency with ~ 2. The qualities of the 
prospective sparticle observations at hadron colliders and linear e^e~ col- 
liders are often very different, with the latters' clean experimental environ- 
ments providing prospects for measurements with better precision. Never- 
theless, Fig. 1121 alreadv restates the clear message that hadron colliders and 
linear e^e^ colliders are largely complementary in the classes of particles that 
they can see, with the former offering good prospects for strongly-interacting 
sparticles such as squarks and gluinos, and the latter excelling for weakly- 
interacting sparticles such as charginos, neutralinos and sleptons. 

Clearly the center of mass energy of any future linear collider is paramount 
towards the supersymmetry discovery potential of the machine. This is seen 
in Fig. El for the benchmark points as more sparticles become observable at 
higher CM energy. We can emphasize this point in general models by plot- 
ting the masses of the two lightest (observable) sparticles in supersymmetric 
models. For example, in Fig. El 36 , a scatter plot of the masses of the 
lightest visible supersymmetric particle (LVSP) and the next-to- lightest visi- 
ble supersymmetric particle (NLVSP) is shown for the CMSSM. Once again, 
points selected satisfy all phenomenological constraints. We do not consider 
the LSP itself to be visible, nor any heavier neutral sparticle that decays 
invisibly inside the detector, such as v ^ vx when v is the next-to-lightest 
sparticle in a neutralino LSP scenario. The LVSP and the NLVSP are the 
lightest sparticles likely to be observable in collider experiments. 

All points shown in Fig. El satisfy the phenomenological constraints dis- 
cussed above. The dark (red) squares represent those points for which the 
relic density is outside the WMAP range, and for which all coloured sparti- 
cles (squarks and gluinos) are heavier than 2 TeV. The CMSSM parameter 
reach at the LHC has been analyzed in [JT]. To within a few percent accu- 
racy, the CMSSM reach contours presented in coincide with the 2-TeV 
contour for the lightest squark (generally the stop) or gluino, so we regard 
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Fig. 12. Summary of the numbers of MSSM particles that may be detectable at var- 
ious accelerators in the updated benchmark scenarios. We see that the capabilities of 
the LHC and of linear e'^e~ colliders are largely complementary. We re-emphasize 
that mass and coupling measurements at e^e~ colliders are usually much cleaner 
and more precise than at hadron-hadron colliders such as the LHC, where, for ex- 
ample, it is not known how to distinguish the light squark flavours. 



the dark (red) points as unobservable at the LHC. Most of these points have 
mNLVSP ^ 1-2 TeV. Conversely, the medium-shaded (green) crosses repre- 
sent points where at least one squark or gluino has a mass less than 2 TeV 
and should be observable at the LHC. The spread of the dark (red) squares 
and medium-shaded (green) crosses, by as much as 500 GeV or more in some 
cases, reflects the maximum mass splitting between the LVSP and the NLVSP 
that is induced in the CMSSM via renormalization effects on the input mass 
parameters. The amount of this spread also reflects our cutoff \Ao\ < 1 TeV, 
which controls the mass splitting of the third generation sfermions. 

The darker (blue) triangles are those points respecting the cosmological 
cold dark matter constraint. Comparing with the regions populated by dark 
(red) squares and medium-shaded (green) crosses, one can see which of these 
models would be detectable at the LHC, according to the criterion in the 
previous paragraph. We see immediately that the dark matter constraint re- 
stricts the LVSP masses to be less than about 1250 GeV and NLVSP masses 
to be less than about 1500 GeV. In most cases, the identity of the LVSP is 
the lighter f. While pair-production of the LVSP would sometimes require a 
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Fig. 13. Scatter plots of the masses of the lightest visible supersymmetric par- 
ticle (LVSP) and the next-to-lightest visible supersymmetric particle (NLVSP) in 
the CMSSM for /x > 0. The darker (blue) triangles satisfy all the laboratory, as- 
trophysical and cosmological constraints. For comparison, the dark (red) squares 
and medium-shaded (green) crosses respect the laboratory constraints, but not those 
imposed by astrophysics and cosmology. In addition, the (green) crosses represent 
models which are expected to be visible at the LHC. The very light (yellow) points are 
those for which direct detection of supersymmetric dark matter might be possible. 

CM energy of about 2.5 TeV, in some cases there is a lower supersymmetric 
threshold due to the associated production of the LSP x with the next light- 
est neutralino X2 jH8j . Examining the masses and identities of the sparticle 
spectrum at these points, we find that Ecu ^ 2.2 TeV would be sufficient 
to see at least one sparticle, as shown in Table 1. Similarly, only a LC with 
EcM > 2.5 TeV would be 'guaranteed' to see two visible sparticles (in addi- 
tion to the X LSP), somewhat lower than the 3.0 TeV one might obtain by 
requiring the pair production of the NLVSP. Points with mj^ygp > 700 GeV 
are predominantly due to rapid annihilation via direct-channel H, A poles, 
while points with 200 GeV < mLvsp ^ 700 GeV are largely due to x-slepton 
coannihilation. 

An EcM = 500 GeV LC would be able to explore the 'bulk' region at 
low (mi/2,"^o), which is represented by the small cluster of points around 
ruLVSP ^ 200 GeV. It should also be noted that there are a few points 
with rriLVSP ~ 100 GeV which are due to rapid annihilation via the light 
Higgs pole. These points all have very large values of mg which relaxes the 
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Table 1. Centre-of-mass energy (in TeV) required to observe one or two sparticles 
at a future LC m the CMSSM and NUHM. 



Model 


sgni/i) 


one sparticle 


two sparticles 


CMSSM 


M > 


2.2 


2.6 




M < 


2.2 


2.5 


NUHM 


M > 


2.4 


2.8 




M < 


2.6 


2.9 



Higgs mass and chargino mass constraints, particularly when mt = 178 GeV. 
A LC with EcM = 1000 GeV would be able to reach some way into the 
coannihilation 'tail', but would not cover all the WMAP-compatible dark 
(blue) triangles. Indeed, about a third of these points are even beyond the 
reach of the LHC in this model. Finally, the light (yellow) filled circles are 
points for which the elastic x-p scattering cross section is larger than 10~* pb. 

Because the LSP as dark matter is present locally, there are many av- 
enues for pursuing dark matter detection. Direct detection techniques rely 
on an ample neutralino-nucleon scattering cross-section. The prospects for 
direct detection for the benchmark points discussed above are shown 
in Fig. 1141 This figure shows rates for the elastic spin-independent and spin 
dependent scattering cross sections of supersymmetric relics on protons. In- 
direct searches for supersymmetric dark matter via the products of annihi- 
lations in the galactic halo or inside the Sun also have prospects in some of 
the benchmark scenarios |39| . 




100 200 300 500 700 1000 100 200 300 500 700 1000 



mj; (GeV) mj; (GeV) 

Fig. 14. Elastic cross sections for (a) spin-independent scattering and (b) spin- 
dependent scattering on protons. Our predictions (blue crosses) are compared with 
those o/ Neutdriver 'lOj (red circles) for neutralino-nucleon scattering. Projected 
sensitivities (a) for CDMS II and CRESST (solid) and GENIUS ^ 
(dashed) and (b) for a 100 kg NAIAD array \44i o.re also shown. 
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In Fig. 1151 we display the allowed ranges of the spin-independent cross 
sections in the NUHM when we sample randomly tan (3 as well as the other 
NUHM parameters . The raggedness of the boundaries of the shaded re- 
gions reflects the finite sample size. The dark shaded regions includes all 
sample points after the constraints discussed above (including the relic den- 
sity constraint) have been applied. In a random sample, one often hits points 
which are are perfectly acceptable at low energy scales but when the param- 
eters are run to high energies approaching the GUT scale, one or several of 
the sparticles mass squared runs negative [JHl- This has been referred to as 
the GUT constraint here. The medium shaded region embodies those points 
after the GUT constraint has been applied. After incorporating all the cuts, 
including that motivated by — we find that the light shaded region where 
the scalar cross section has the range 10~^ pb > asi ^ lO""'^" pb, with some- 
what larger (smaller) values being possible in exceptional cases. If the — 2 
cut is removed, the upper limits on the cross sections are unchanged, but 
much lower values become possible: asi <C 10~^^ pb. The effect of the GUT 
constraint on more general supersymmetric models was discussed in |47| . 

The results from this analysis 131 for the scattering cross section in the 
NUHM (which by definition includes all CMSSM results) are compared with 
the previous CDMS and Edelweiss H^I bounds as well as the recent 
CDMSII results [SUI in Fig.^] While previous experimental sensitivities were 
not strong enough to probe predictions of the NUHM, the current CDMSII 
bound has begun to exclude realistic models and it is expected that these 
bounds improve by a factor of about 20. 

This work was partially supported by DOE grant DE-FG02-94ER-40823. 
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